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We present the first four-flavour lattice calculation of the leading-order hadronic vacuum- 
7— I polarisation contribution to the anomalous magnetic moment of the muon, a^ p , and the hadronic 

running of the QED coupling constant, Aa^ D (Q 2 ). In the heavy sector a mixed-action setup is 
y— I employed. The bare quark masses are determined from matching the K- and D-meson masses to 

their physical values. Several light quark masses are used in order to yield a controlled extrapola- 
tion to the physical pion mass by utilising a recently proposed improved method. We demonstrate 
that this method also works in the four-flavour case. 
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1. Introduction 

The anomalous magnetic moment of the muon, , serves as a benchmark test of the standard 
model (SM). It has been measured very accurately [1, 2] and can be computed precisely within 
the SM. A comparison between the experimental result for and the SM prediction reveals a 
discrepancy of more than three standard deviations which has persisted for many years now and 
has been confirmed by computations of a number of groups, see, for example, the review [3]. The 
question is whether this discrepancy originates from some effect missing in the experimental or 
theoretical determination of or whether it points to physics beyond the SM. 

A key ingredient in the calculation of is the leading-order hadronic vacuum-polarisation 
contribution, a^ v , which presently is the largest source of uncertainty in the theoretical compu- 
tation of a^, since the QED and electroweak contributions have been computed very accurately 
employing perturbation theory, see [4, 3] and references therein. As a^ v is intrinsically nonpertur- 
bative, a lattice QCD computation of this observable is highly desirable. The currently accepted 
SM values for this quantity are obtained mainly from the investigation of e + e~ scattering and z 
decay data. In a recent study [5, 6], using two flavours of mass-degenerate quarks, a modified 
method to compute a|7 P has been introduced resulting in a determination of a h ^ iaht with a precision 
of a few percent. 

Besides a^ p , the leading-order QCD contribution to the running of the QED coupling constant, 
AoCqed' also requires the hadronic vacuum-polarisation function and can likewise be investigated 
once this function is known. As an important input parameter to SM calculations, the QED coupling 
constant needs to be known very precisely in order to facilitate high precision tests of the SM at 
any future linear collider [7]. 

Below we report preliminary results of the first lattice calculation of a^ v and Aagg P D with four 
quark flavours. The contributions of the top and bottom quarks are negligible at the current level of 
accuracy. So having four flavours allows for an unambiguous comparison to the dispersive analysis 
of a|7 P an d direct use in forming the SM prediction for itself. The inclusion of the charm quark 
is essential because its contribution is of the order of aj7 P ' charm ^ 100 x 10~ n [9], which is larger 
than the currently quoted uncertainty of the difference between the experimental and the SM results. 
Furthermore, the order of magnitude of the charm quark contribution to is the same as that of 
the hadronic light-by-light contribution [10]. Thus lattice calculations with a dynamical charm 
quark are necessary for computing with a precision that matches the experimental accuracy. 

The computation of a^ v and Ao!q£ P d follows closely the strategy of refs. [5, 6] using improved 
lattice definitions of these quantities. We demonstrate in this work that this new method continues 
to work well even for the four-flavour case and again results in a mild quark mass dependence 
leading to an accurate determination of a^ p and Ao!q£ P d . 

Our calculations are based on the configurations with four dynamical quark flavours generated 
by the European Twisted Mass Collaboration (ETMC) [11, 12]. These sets of configurations are 
obtained at different values of the lattice spacing and several lattice volumes, thus enabling us to 
estimate discretisation and finite size effects as systematic uncertainties of our lattice calculation. 
In addition, at each value of the lattice spacing, configurations exist at several values of the pion 
mass, ranging from 230 Me V < m % < 450 Me V. 
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2. The lattice calculation 



The leading-order hadronic contribution to the muon's anomalous magnetic moment, aj7 P , can 
be computed directly in Euclidean space-time [13] 

The renormalised vacuum-polarisation function is given by I1r(<2 2 ) = Tl(Q 2 ) — 11(0). n(<2 2 ) can 
be obtained from the hadronic vacuum-polarisation tensor knowing its Lorentz structure 

iW<2) = / d 4 xe^-y\jfr( x )J7(y)) = (QnQv - Q 2 s,vMq 2 ) , (i.i) 



where 



2_ 1_ 2_ 1_ 

J T(x) = ^u{x)y lJL u{x) - -d{x)y^d{x) + -c(x)y^c(x) - -s(x)Yns(x) (2.3) 



is the electromagnetic vector current. 

On the lattice we use the conserved (point-split) vector current 

= ^(^+A)(i+r^WGeizW 

- ^(x)(l-7 M )t/ M (x)e el ^(x+A)) • (2.4) 

X(x) denotes a fermion doublet, either the light or the heavy one, in the so-called twisted basis 
and <2ei = diag(|, — |) denotes the electric charge matrix. Resorting to Osterwalder-Seiler valence 
quarks [14, 15] in the heavy sector admits a straightforward construction of the conserved currents 
also for the strange and the charm quark. This is beneficial since in this way we can rely on the 
vector Ward-Takahashi identity for all contributions and also avoid renormalisation. 

In order to have a smooth function to perform the integral in eq. (2. 1), the vacuum-polarisation 
data obtained at discrete lattice momenta is fit for each flavour to the following functional form 

2 

U(Q 2 )=g 2 J M v +b + b l Q 2 . (2.5) 
Q l + mif 

The first term is the contribution of a narrow- width vector meson with mass my and electromagnetic 
coupling gy, which are determined directly in our calculation. The remaining terms parametrise 
any deviations from this form. The results reported below are obtained using uncorrelated fits to 
determine gv and my from the zero-momentum current correlators. In these proceedings we do not 
perform alternative fits such as the Pade approximants suggested in [17]. We also do not provide 
an estimate of systematic effects by allowing additional terms in eq. (2.5). These aspects will be 
addressed in a forthcoming publication. 

Once I1r(<2 2 ) is known, it is straightforward to compute the leading-order hadronic contribu- 
tion [7] 

A «Q v E p D (e 2 ) = 4nabii R (e 2 ) (2.6) 

which influences the running of the fine structure constant according to [16] 

^ rfl ° i-AMW <2 " 7) 
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Here, Ofo is the value at Q 2 = 0, <Xq 

For the lattice calculation of aj7 P discussed here, we use the modified definition from [5] 



hvp 



«0 



Q 2 



-w 



e 2 # P V 



H 2 



1117 



n R (e 2 



(2.8) 



which goes to a^ v for the light pseudoscalar mass mps assuming its physical value, i.e. m n , since in 
this case also H — > // p h ys , with the choice of the hadronic scale, H, discussed below. Analogously, 
we determine Ao!q£ P d (<2 2 ) from [6] 



A«Q V E P D (e 2 )=47ra n R 2 : 



H 1 



(2.9) 



There are several possible choices for H. Below we give results for H = // p h ys = 1 (the standard 
choice) and H = rav('Wps), # P hys = »ip (improved choice). Here mv(n?ps) is the light vector meson 
mass as measured on the lattice at unphysical pseudoscalar masses nips whereas m p is the physical 
value of the p-meson mass. See refs. [5, 6] for a more detailed discussion. 
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Figure 1: Light quark contribution to a^ p as a function of the squared light pseudoscalar mass from our 
Nf = 2 + 1 + 1 computations. The lower set of data points corresponds to using the standard definition, H = 
^phys = 1 in eq. (2.8). The upper data points are obtained using the improved method, setting H = m\r(mps) 
and Z/phys = fn p . The open square represents the two-flavour result of refs. [5, 6]. The dark grey error band 
corresponds to the quadratic fit (solid green line) and the light grey one belongs to the linear fit (dashed black 
line). 

The light quark contribution to a^ p is depicted in fig. 1 . In contrast to the standard choice the 
improved lattice definition of a^ v shows a weaker pion mass dependence. Comparing the previous 
calculation with only two dynamical flavours of light, mass-degenerate quarks [5] to the one having 
four flavours in the sea quark sector, the results for ajf^ extrapolated to the physical point are 
found to be in full agreement for both setups. This demonstrates that the effects of a dynamical 
strange and charm quark on the light quark contribution to a]j vp are small, as expected. 

In fig. 2 we show the effect of adding the strange quark contribution and compare to the pub- 
lished results of other groups, [18, 19]. The figure demonstrates that when using the standard 
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definition of a^ p on the lattice, all groups agree reasonably well. However, when using the im- 
proved definition, the pion mass dependence of a^ p is much flatter allowing for a better control of 
the extrapolation to the physical pion mass. 
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Figure 2: Comparison of the results for aj? p having the light up and down quarks as well as the sttange 
quark in the valence sector for different collaborations. Only the upper set of data points for twisted mass 
fermions are obtained from the modified definition, H = my(»ips) and // p h ys = mp in eq. (2.8). Data are 
from [19] (clover-improved) and [18] (domain wall). The open square represents the standard model value 
obtained from the dispersive analysis of ref. [8]. Concerning the error bands of the fits, the same comments 
as in fig. 1 apply. 

Finally, we show in fig. 3 the full four-flavour contribution to a^ p . Comparing with fig. 2, 
we see that including the charm quark indeed leads to a contribution of the expected order of 
magnitude. We perform both a linear and a quadratic extrapolation to the physical point. Within 
just the statistical errors, we find reasonable agreement with the dispersive result [8], as shown 
in fig. 3. At the current precision, we do not observe any statistically significant finite-size or 
lattice discretisation effects, however, the impact of disconnected contributions or unitary-violating 
effects due to the mixed-action setup in the heavy sector have not been examined yet. Furthermore, 
the effect from the p -meson not being able to decay to two pions for all but one ensemble has 
to be investigated. The details of the final extrapolation to the physical point and the associated 
systematic uncertainties will be presented elsewhere. 

Fig. 4 shows our results for AotQg P D at a typical hadronic scale of Q 2 = 1 GeV 2 obtained from 
the same ensembles used to determine a^ p . Similar to a^ p , we find reasonable agreement with the 
dispersive result [7] to within the statistical uncertainties. A full determination of the systematic 
uncertainty for the extrapolated value will also be presented in a later publication. 

3. Conclusion 



We have pointed out that the charm quark contribution is necessary for achieving a lattice 
computation of a^ p with a precision that is comparable to the experimental one. Furthermore, 
we have shown that also in the case of a four-flavour calculation of the leading-order hadronic 
contribution to the muon anomalous magnetic moment and to the running of the electromagnetic 
coupling constant, the improved method of ref. [5, 6] continues to work well. We have so far not 



5 



aj7 P an d A&QEY)f>~om Nf = 2 + 1 + 1 twisted mass fermions Grit Hotzel 



7e-08 
6e-08 
5e-08 
i 4c-08 
3c-08 - 

CO 

2e-08 - 
lc-08 - 




(I 



1 ! 



Preliminary 



dispersive analysis 
a = 0.078 fm, L = 1.91m 
a = 0.078 fm, L = 2.5 fm 
a = 0.061 fm, L = 1.9 fm 
a = 0.061 fm, L = 2.9 fm 



0.05 



0.1 



0.15 

teeV 2 l 



0.2 



0.25 



Figure 3: Full four-flavour contribution to a^ vp as a function of the squared light pseudoscalar mass from 
our Nf = 2 + 1 + 1 computations. Concerning the upper and lower data sets and the error bands of the fits, 
the same comments as in fig. 1 apply. The open square represents the standard model value obtained from 
the dispersive analysis of ref. [8]. Note that in this case of four flavours there is no ambiguity in determining 
the contribution to a^ v from the dispersive analysis. 
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Figure 4: Full four-flavour contribution to Ao!Qgp(l GeV 2 ) as a function of the squared light pseudoscalar 
mass from our Nf = 2 + 1 + 1 computations. The lower set of data points corresponds to the use of the 
standard definition, H = H p h ys = 1 in eq. (2.9). The upper data points are obtained using the improved 
method, setting H = m\(mps) and // p h ys = '«p- Concerning the error bands of the fits, the same comments 
as in fig. 1 apply. The open square represents the standard model value obtained from the dispersive analysis 
of ref. [7]. Note that in this case of four flavours there is no ambiguity in determining the contribution to 
AaQg P D (l GeV 2 ) from the dispersive analysis. 



performed a comprehensive analysis of the systematic uncertainties. In particular, the disconnected 
contributions originating from the strange quark might be non-negligible. Additionally, it might 
be necessary to take isospin breaking effects into account as the precision of our computation 
improves. 

We conclude by emphasising that the strategy followed here can also be applied to other ob- 
servables such as the Adler function, the corrections to the energy levels of muonic and ordinary 
hydrogen, and the weak mixing angle, as discussed in [6]. 
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